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Abstract 

Assuming that neutrinos acquire radiative seesaw Majorana masses through their 
interactions with dark matter, i.e. scotogenic from the Greek 'scotos' meaning dark- 
ness, and using the non-Abelian discrete symmetry A4, we propose a model of neu- 
trino masses and mixing with nonzero #13 and necessarily large leptonic CP violation, 
allowing both the normal and inverted hierarchies of neutrino masses, as well as quasi- 
degenerate solutions. 



In 2006, a one-loop mechanism was introduced [T] linking neutrino mass with dark matter. 
The idea is very simple. The standard model of particle interactions is extended to include 
a second scalar doublet (i] + ,r]~) which is odd under an exactly conserved Z 2 symmetry [2], 
as well as three neutral fermion singlets iVj which are also odd under Z 2 . This requirement 
immediately allows the possibility of having the real (or imaginary) part of rf as a dark- 
matter candidate, which was first pointed out also in Ref. pQ. As shown in Fig. 1, this 
results in the radiative generation of seesaw Majorana neutrino masses from dark matter, 
i.e. scotogenic from the Greek 'scotos' meaning darkness. 
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Figure 1: One-loop generation of scotogenic Majorana neutrino mass. 

The non-Abelian discrete symmetry v4 4 was introduced [SI SI E] to achieve the seemingly 
impossible, i.e. the existence of a lepton family symmetry consistent with the three very 
different charged- lepton masses m e , m M , m T . It was subsequently shown [6J to be a natural 
theoretical framework for neutrino tribimaximal mixing, i.e. sin 2 #23 = 1, tan 2 9 12 = 0.5, and 
#13 = 0. This pattern was consistent with experimental data until recently, when the Daya 
Bay Collaboration reported [7] the first precise measurement of 9 13 , i.e. 

sin 2 29 13 = 0.092 ± 0.016(stat) ± 0.005(syst), (1) 

followed shortly [H] by the RENO Collaboration, i.e. 

sin 2 29 13 = 0.113 ± 0.013(stat) ± 0.019(syst). (2) 
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This means that tribimaximal mixing is not a good description, and more importantly, 
leptonic CP violation is now possible because #13 7^ 0, just as hadronic CP violation in the 
quark sector is possible because V u b 7^ 0. 

Recently, it was shown [9] that A4 is still a good symmetry for understanding this pattern, 
using a new simple variation of the original idea [6]. In that proposal, neutrinos acquire 
Majorana masses through their direct interactions with Higgs triplets. We study here instead 
the corresponding scenario with the radiative mechanism of Fig. 1. 

The symmetry A4 is that of the even permutation of four objects. It has twelve elements 
and is the smallest group which admits an irreducible three-dimensional representation. Its 
character table is given below. The basic multiplication rule of A4 is 

3 x 3 = 1 + 1' + 1" + 3 + 3. (3) 

As first shown in Ref. [3], for (z/j, Zj) ~ 3, l\ ~ 1,1', 1", and $j = (0°, ~ 3, the charged- 
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Table 1: Character table of A4. 



lepton mass matrix is given by 
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/ a/3, we then obtain 
(5) 
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where m e = fiv, m M = /2U, m T = f$v. The original A4 symmetry is now broken to the 
residual symmetry Z3, i.e. lepton flavor triality [ID] , with e ~ 1, fi ~ w 2 , r ~ u. This is a 
good symmetry of the Lagrangian as long as neutrino masses are zero. Exotic scalar decays 
are predicted and may be observable at the Large Hadron Collider (LHC) in some regions 
of parameter space [HI [12] . 

To obtain nonzero neutrino masses, we assign rj ~ 1 and Ni ~ 3 under A4. We also add 
the scalar singlets Oi ~ 3 with nonzero (o"j). The resulting 3x3 Majorana mass matrix for 
Nj is then 



M 



N 



I A F E\ 

FAD 
K EDAj 



(6) 



which is the analog of 



fa f e\ 

M u = f a d , (7) 
\ e d a J 

considered in Ref. [9]. (A better way to enforce Eq. (6) is to postulate gauged B — L and 
assume complex neutral scalars which transform as 1, 3 under A4, in complete analogy to 
the scalar triplets of Ref. [9].) Instead of enforcing E = F = which is required for tribi- 
maximal mixing, we assume here that F = —E which may be maintained by an interchange 
symmetry P Q2] . 



Consider now the tribimaximal basis, i.e. 



fu e \ ( J2/3 l/\/3 



-1/VQ l/y/3 -1/V2 

V-1/V6 l/y/3 1/V2 J 



Since are connected to ^1^,3 through the identity matrix, we find 



/A + D \ 

A C 
K C A-Dj 



(9) 



where C = (E - F)/y/2 = V2E. 
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The diagram of Fig. 1 is exactly calculable from the exchange of Re(r]°) and lm(j]°) and 
is given by [1] 

,2 ™2 



/ a a \ \ - hikhjkMk 



m\ m 2 R in 
m 2 R - Ml n Ml~mj-M 2 k ^Ml 



1 InS 



(10) 



where J2khik(hjk)* = \h\ 2 5ij, and m^/ are the masses of y/2Re(i]°) and v^2bm(?7 ), respec- 
tively. In the limit m 2 R — m 2 = 2X^v 2 is small compared to m\ = {m 2 R + mj)/2, and 
ml « Mf , Eq. (10) reduces to 
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In the tribimaximal basis of Eq. (9), we then have 
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(13) 



The neutrino mixing matrix U has 4 parameters: S12, S23, S13 an< l ^cp [14J . We choose 
the convention U T i,U T 2,U e 3,U^ — > —U T i,—U T 2,—U e s,—U^ to conform with that of the 
tribimaximal mixing matrix of Eq. (8), then 



where m' 12 3 are the physical neutrino masses, with 
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(14) 



m' 2 + + Am^ (normal hierarchy), 



vn' 2 + A777J1/2 — Am 2 2 (inverted hierarchy). 



(15) 
(16) 
(17) 
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We now diagonalize M^' 2 ' 3 ^ using 

/e <a imi \ 

U e M { J' 2 ' 3) Uj = e ia 2m' 2 , (18) 

v e ia{ im' 3 J 

from which we obtain U' = UtbUJ '. To obtain U with the usual convention, we rotate the 
phases of the fx and r rows so that U^ 3 e~ ia ^ 2 is real and negative, and U! r3 e~ ta ^ 2 is real 
and positive. These phases are absorbed by the \x and r leptons and are unobservable. We 
then rotate the i/ lj2 columns so that U' el e~ ia ^ 2 = U el e ia "l 2 and U' e2 e~ ia ^ 2 = U e2 e ia ^ 2 , where 
U e i and U e 2 are real and positive. The physical relative Majorana phases of v 12 are then 
<*i,2 — Qj'i 2 + a i,2- The three angles and the Dirac phase are extracted according to 

tan 2 # 12 = \U' e2 /U' el \ 2 , tan 2 ^ 3 = \U^/U' tZ \\ sin^e"^ = C^e"^/ 2 . (19) 

The effective Majorana neutrino mass in neutrinoless double beta decay is then given by 

m ee = \U 2 1 e iai m' 1 + U 2 2 e ia2 m 2 + U 2 3 m' 3 \. (20) 



In Eq. (9), let A be real and positive by convention, then both C and D may be complex, 
i.e. C = Cr + iC[ and D = D R + iDj. The 2x2 matrix of Eq. (13) can be solved exactly 
to yield 

CrDj - CjDr 



tan0 
tan 29 



C R {2A-D R )-C I D I ' 

2[AA 2 C 2 R - AACr{CrDr + CjDj) + (C 2 R + Cf)(D 2 R + D 2 )} 1 ' 2 
2AD R - (D 2 R + D 2 ) 



with 



e ia 2M 2 = cos 2 9 A + 2 sin 9 cos Be* C + sin 2 9e 2i<t> (A- D), 
e ia 'sM 3 = cos 2 9{A-D) - 2 sin 9 cos 9e^C + sin 2 9e~ 2l(t> A. 



The corresponding U' elements are 
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If we absorb the scale factor \ 5 h 2 v 2 /87r 2 into the parameters A, C, D as well as m , then 
the mass eigenvalues of Eq. (11) are given by 



m ' k = w k 



In— £-- 1 
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(27) 



which are the ones used in Eqs. (14) and (18). Since m is an unknown, having to do with 
the dark-matter scalar mass, we fix it by requiring Mi /mo = 10, where M\ = \A + D\. If we 
input the five parameters A, Cr, Cj, Dr, Dj, we will obtain m\ 2 3 as well as the three mixing 
angles and the three CP phases. For our numerical analysis, we set 

Am 2 21 = 7.59 x 1(T 5 eV 2 , Am 2 2 = 2.45 x 1(T 3 eV 2 , (28) 

and vary # 13 in the range 

sin 2 20i3 = 0.05 to 0.15. (29) 

Following Ref. [9], we look for solutions with sin 2 2# 2 3 = 0.92 and 0.96. Whereas only normal 
hierarchy is allowed in the model of Ref. [9] , we find solutions for both normal and inverted 
hierarchies, as well as quasi-degenerate solutions, as detailed below. 

The predictions of this model regarding mixing angles are basically the same as in Ref. [9] 
for the special case of b = there. Using Eqs. (19), (25), and (26), we find 

tan 2 ^ = 1 ~ 3sin2gl3 , (30) 



y^2 sin 0i3 cos <j> \ . 2 sin 2 613 sin 
V /l - 3& ' in2e l3/ l-3sin 2 9 13 



2 " 2. 



tan 6>23 = y -2 • (31) 

/ -I 1 y/2 sin 613 cos <f> \ , 2 sin 2 ^13 sin 2 
V V!-3 8in 2 0i3/ l-3sin^ 13 

The conventionally defined Dirac CP phase is given by Sep = <f) + ot' 3 /2, where a' 3 is defined 
in Eq. (18) and depends on the specific values of Eq. (9). For sin#i 3 = 0.16, corresponding 
to sin 2 2013 = 0.1, this predicts tan 2 9 12 = 0.46. If Im{C) = 0, then 5cp = a' 3 = 0, so this 
would predict sin 2 26*23 = 0.80 which is of course ruled out. Using sin 2 26*23 > 0.92, we find 
in this case I tan</>| > 1.2. 



For each of the two values sin 2 26*23 = 0.92 and 0.96, we obtain 5 representative solutions, 
all as functions of sin 2 26 ) i3. Using Eq. (30), we plot sin 2 29±2 versus sin 2 26*13 in Fig. |2j The 
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Figure 2: sin 2 26*12 versus sin 2 2# 13 . 



characteristic features of the 5 solutions are listed in Table 2. For Im(D) = 0, we find one 
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Table 2: Five representative solutions. Three have Im(D) = 0, and two have Im(D) = 
Re(D). NH denotes normal hierarchy of neutrino masses, IH inverted, and QD quasi- 
degenerate. The values of | tan<5cp| and m ee (in eV) are for sin 2 26*23 = 0.96 and sin 2 26*i3 = 
0.10. 

solution for inverted ordering of neutrino masses, and two solutions for normal ordering (one 
of which is quasi-degenerate). For Im(D) = Re(D), we again find one solution for inverted 
ordering, but the only solution for normal ordering is quasi-degenerate. 

In Fig. 3 we show the physical neutrino masses 2 3 and the effective mass in neutrinoless 
double beta decay m ee (in eV) as well as the model parameters (in eV -1 ) for solution (I) in the 
case sin 2 26*23 = 0.96. In Figs. 4-7 we show the same quantities for solutions (II), (III), (IV), (V) 
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in the cases of sin 2 2#23 = 0.92, 0.96, 0.92, 0.96 respectively. Finally we show in Fig. 8 the 
values of |tan<5cp| for all 5 solutions in the case of sin 2 2^23 = 0.92. It is claer that at 
sin 2 26*13 = 0.10, large | tax\.5 C p\ is predicted. 
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Figure 3: A 4 parameters and the physical neutrino masses and effective neutrino mass m ee in 
neutrinoless double beta decay for the inverted hierarchy with Im(D)=0 and sin 2 2#23 = 0.96. 




Figure 4: A4 parameters and the physical neutrino masses and effective neutrino mass m ee in 
neutrinoless double beta decay for the inverted hierarchy with Im(D)=Re(D) and sin 2 2#23 = 
0.92. 
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Figure 5: A 4 parameters and the physical neutrino masses and effective neutrino mass m ee in 
neutrinoless double beta decay for the normal hierarchy with Im(D)=0 and sin 2 2^23 = 0.96. 




Figure 6: parameters and the physical neutrino masses and effective neutrino mass m ee 
in neutrinoless double beta decay for quasi-degenerate neutrino masses with Im(D)=0 and 
sin 2 2003 = 0.96. 
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Figure 7: A4 parameters and the physical neutrino masses and effective neutrino mass m ee 
in neutrinoless double beta decay for quasi-degenerate neutrino masses with Im(D)=Re(D) 
and sin 2 26 23 = 0.96. 




